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5} 1 Abstract 

Q_i. Modifications of the mass spectrum the quarkonium induced by its vir- 

tual dissociation into a pair of heavy mesons is considered. Coupling be- 
CN \ tween quark and mesonic channels results in noticeable corrections to spin- 

dependent mass splitting. In particular, the observable hierarchy of mass 
splittings in the Xc Xb and x'b multiplets is reproduced. 
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Quarkonium is the testing ground for various theoretical approaches and, 
especially, of potential models. In particular, analysis of the quarkonium 



(N 



CO 



mass spectrum is used for investigating properties of the confining potential. 
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Phenomenological potentials were reasonably successful in describing gross 
features of the quarkonia spectra: the level spacing and ordering . More sub- 
"r* ■ tie details of the mass spectra, in particular, spin- dependent mass splitting 

S^h! still remains a challenging problem. Time-honored approach to the analysis 

of spin-dependent effects in quarkonia consists in introducing an interquark 
QED-type potential but with more generic vertex structure pj: 



V = TtViq 2 )^, (1) 

where V(q 2 ) is the propagator of the exchanged particle and Tj is the coupling 
vertex (equal e.g. to 1 or 7^ in the case of scalar and vector-type coupling 
respectively) . In the nonrelativistic limit the expansion of (0) in the inverse 
powers of the heavy-quark mass yields the sum of central and spin-dependent 
components 



V = V C + 



(Si + S 2 ) L] V lo + Sxz ■ V T + (Sx ■ S 2 ) V ss (2) 



with Vlo, Vt, Vss standing for the spin-orbit, tensor and spin-spin potentials 
respectively. In eq. (H) Sip are the quarks' spins, L is their orbital momentum 
and 5*12 = 12(5*1 • r )(S 2 ■ r) — 4(Si • 52). The spin-structure of the interquark 
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potential (0) can be translated in the mass splitting of the triplet P states 
as follows: 

m( 3 P ) = m — 2a — 4b 
m( 3 Pi) = m — a + 2b 
m( 3 P 2 ) = m + a — 26/5 

where fh is the c.o.g. mass of the triplet (weighted with 2 J + 1 factor), 
a and b are the averaged values of the Vlo and Vr potential components 
respectively. The Vss component of the potential controls mass splitting 
between the triplet and singlet P states. The combined action of spin-orbit 
and tensor forces is parameterized in the form of the ratio 

m( 3 P 2 )-m( 3 P) 

m( 3 Pi) - m( 3 P ) 1 ' 

The masses of the 3 Pj charmonium states as determined by Crystal Ball 
yield 

R( Xc ) = 0.478 ± 0.01 (4) 

with corresponding ratios for l 3 Pj and 2 3 Pj states of the 66-quarkonium 
reading respectively 



R( Xb ) = 0.664 ± 0.038 (5) 



and 



R( x ' b ) = 0.576 ± 0.014 (6) 
Note the hierarchy of P's as given by the measured quarkonia masses: 

R e * p (x b )>R exp (x' b )>R e * p (xc) (7) 

At the same time potential models involving the sum of the the scalar (con- 
fining) and vector (residual one-gluon exchange) potentials (c/(|) ) yielded 
controversial results predicting R(xb) > R(Xc) and R(x'b) slightly larger or, 
at best, approximately equal to R(xb) i n the qualitative disagreement with 
(0) (see § for the details): 

R(x' b )>R(Xb)> R(Xc) (8) 

Observed disagreement hints at the existence of additional effects. A solution 
to the problem was sought by introducing some unconventional interactions, 
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in particular, a pseudoscalar component of the quark-antiquark potential 
with T = 7 5 in ([!]) [|J. Another issue might be an account of relativistic 
effects which are more pronounced in cc as compared to their bb counter- 
parts. Currently no convincing solution to the problem in the framework of 
the potential-based approach seems to exist. The controversy casts doubts 
on the universally adopted conclusion about scalar Lorentz-transformation 
properties of the confining potential. 

In the present paper we consider an alternative mechanism affecting mass- 
splitting in the quarkonia systems which might give an insight to the solu- 
tion of the challenging problem. It is coupling of quarkonia states to pairs 
of heavy quark-mesons HH. Investigation of the effects emerging from such 
coupling was pioneered by 0, where a unified description of the heavy quark- 
antiquark binding and of the quarkonium coupling to heavy mesons was 
given. However, the calculations made based heavily on the specific form of 
the confining potential. The latter was assumed to have Lorentz- vector trans- 
formation properties in disagreement with the current experimental data on 
both the spectroscopy of cc states and on their decay amplitudes . The lat- 
ter circumstance substantiates revisiting effects which emerge from coupling 
between quarkonia states and mesonic channels. 

Coupling of quarkonium (QQ) bound states to heavy-meson pairs (HH) 
results in the mass shift of the former because of the virtual dissociation 

QQ - (Qq) + {Qq) - QQ (9) 

The mass shift emerges as a matter of fact due to an interaction between 
heavy QQ and (virtual) light qq quark pairs. The situation, however, is 
somewhat tricky. Indeed, in the adiabatic approximation the presence of 
the light-quark pairs does not affect (at least) the long-range part of the 
potential operating between heavy quark and antiquark. More precisely, 
its linear form is maintained, and only the tension is altered. The latter, 
however, is a phenomenological parameter determined from the experimental 
data implying that in the uiq — > oo limit coupling to mesonic channels 
does not result in an observable effect J7J. This conclusions refers to the 
spectrum gross features controlled by the string tension. At the same time 
spin-dependent mass splitting is scaled by a 1/mg factor, and non- vanishing 
effect are possible. 

Consider in more detail the mechanism which produces effective mass 
splitting of the Pj triplet. The lightest Qq states (and their antiparticles) 
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materialize in the form of pseudoscalar and vector mesons (denoted hereafter 
H and H* respectively). Be H and H* mass-degenerate, coupling between 
QQ and qq would not result in mass splitting of quarkonium states. However, 
because of the spin Q — q interaction the H* vector meson proves to be heavier 
than its pseudoscalar H counterpart. It implies that the mass shift of the 
quarkonium states is scaled by the width of a HH continuous spectrum band 
up to the effective mass >(m#* — wig). The contribution coming from HH 
pairs with masses exceeding tuh + inn* is canceled by similar contribution 
of HH* and H*H* heavy-meson states. 

A quantum-mechanical consideration of the problem allows making some 
qualitative conclusions concerning the expected effect. First, quarkonium 
states, depending on their spin-parity quantum numbers, couple to different 
HH states whereas some of them do not have a HH match at all. Indeed, the 
3 Po an d 3 -P2 quarkonium states couple to the HH pairs with J n = + and 2 + 
respectively, while the 3 Pi level is uncoupled from the (pseudoscalar) mesonic 
sector. The 1 P\ state does not couple to HH mesons either. Selective cou- 
pling results in different mass shifts of the QQ states, and, eventually, in their 
mass splitting. Second, the HH continuous spectrum lies above the quarko- 
nium bound states. Spin-independent potentials reproduce the gross features 
of the spectrum accurately enough to assume that the spin-dependent forces 
may be treated as a perturbation. The second-order perturbation in energy 
[§] equals 

\y I 2 

£( 2 ) = V (10) 

*n 2^ „(0) p (0) ^ 

where V nm is a perturbing interaction and E^ is the (unperturbed) energy 
of the n-th state. Eq . ([10|) shows that the mass of the lowest state E$ (i.e. 
the quarkonium) decreases because of channel coupling (since Eq < E$ ). 

Combining these observations we can conclude that the masses of both 
the 3 Po an d 3 ^2 states will diminish because of their coupling to mesonic HH 
states, whereas the mass of the 3 Pi state will remain unaffected. Inspection of 
(0) shows immediately that the considered mechanism results in the decrease 
of the R ratio. 

Next observation is that the magnitude of the effect depends on the width 
of the energy gap separating quarkonium states from the continuous HH 
spectrum: the smaller is the gap, the more pronounced will be the mass 
shift. We infer then that the R(xb) ratio (l 3 Pj bottomonium states) will 
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be affected less than the R(x'b) ratio (2 3 Pj bottomonium states). Compar- 
ing these qualitative predictions to the experimental data (^) we conclude 
that the considered mechanism (coupling between quarkonium and mesonic 
states) gives a welcome trend to the modified theoretical results. 

We proceed now to quantitative estimates of the effect. To this end 
a dynamical model describing the coupling mechanism under consideration 
is to developed. We treat quarkonium (QQ) and mesonic (HH) states as 
the components of a two-channel system labeling them channels 1 and 2 
respectively. Then the Schrodinger equations describing the dynamics of the 
system for each J 71 " state can be written as follows 

Tx^i + V 11 <pi+ V 12 <f)2= #i0i ni N 

Here T is the kinetic energy in the corresponding channels and V is the 
interaction potential. The Vn component is a (linear + Coulombic) po- 
tential operating between heavy quarks. Interaction between if-mesons is 
neglected: near-threshold phenomena under consideration involve small rel- 
ative momenta and, hence, large inter-particle distances. At the same time, 
the long-range interaction - 7r— meson exchange - between two pseudoscalar 
particles is absent. Anyway, this approximation does not alter the qualita- 
tive conclusions made above. Finally, E\ and E 2 are the channel energies, 
differing by the mass gap in the corresponding state: E 2 — E\ = 2m H — rnqq. 
The non-diagonal component of the potential V\ 2 (= V 2 \) describes the dy- 
namics of the quarkonium dissociation (QQ) — > HH, or, stated differently, 
the dynamics of the qq production in the field of heavy quarks. Forbidding 
complexity of the problem necessitates introducing some simplifying approx- 
imations. 

Given the large mass difference between heavy and light quarks, we make 
use of the adiabatic approximation. The latter molds in the following as- 
sumptions. First, the position of the heavy quarks is not affected by the 
process of the qq production. Second, the combined spin of heavy quarks in 
a pair of heavy mesons is equal to zero. It implies that in all of the tran- 
sitions considered the spin of the QQ pair is subject to one and the same 
change from S = 1 to S = 0. Finally, c.o.g. of a heavy meson coincides with 
good accuracy with the position of the heavy meson residing in it. All the 
assumptions made above can be expressed concisely in the following form of 
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the non- diagonal component: 



V 12 = (3 5(r-r') 



(12) 



where (3 is a coupling strength constant, f and r ' are the separation between 
heavy quarks prior and after the qq production (i.e. within the quarkonium 
state and in the HH pair respectively). The light-quark spin operator and 
coordinate variables are assumed to be present and are dropped for notational 
brevity. 

In the framework of the model the strength constant (3 can be determined 
without going into details of the complicated mechanism of the light-quark 
pair production. The same coupling between channels governs decays of the 
quarkonium states lying above the HH threshold, i.e. -^(3770) and T(10580) 
in the case of charmonium and bottomonium respectively. They are known 
to decay predominantly into the DD and BB pairs with the widths equal 
to 83.9±2.4 and 21±4 MeV respectively. The quantum numbers of the bot- 
tomonium state are known (J 77 = 1~) translating into the 4 3 Si state of the 
bb pair. The t/>(3770) quantum numbers are not determined; following the 
potential model prescriptions we adopt that it is the l 3 Di state of the cc pair 
[0- |]- Using the channel-coupling potential (|T2|) we calculate the value of 
the (3 strength coupling constant. It proves to be equal to 153 MeV and 114 
MeV for the charmonium and bottomonium systems respectively. The agree- 
ment between two values of (3 is reasonably good signaling the validity of the 
approximations made. Indeed, the paradigm of the non-perturbative QCD 
assumes the light-quark pair production as a result of the string break-up. 
Flavor-blindness of the confinement forces (materialized as a string) prompts 
then that in the infinite heavy-quark mass limit the values of j3 for two types 
of the quarkonium should the same. Given possible 1/rriQ corrections, the (3 
values may be considered as coinciding in line with the theoretical expecta- 
tion. 

Having determined the potential V\2 connecting the quarkonium and the 
mesonic channels we proceed to calculating the matrix element of the transi- 
tion between them. The matrix element controls the value of the mass shift 
(see ([Top). Within the approximations made it reduces to the overlap integral 
of the quarkonium wave function ip®® and of the wave function ip^ H of two 
non-interacting heavy mesons: 




(13) 
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The subscript v labels the states of the HH continuous spectrum, to this end 
we use the c.m.s. momentum of heavy mesons iy = k). Then the expression 
for the quarkonium-state mass shift ([II]) can be cast in the form 



Am QQ 



m QQ ~ m HH 



dk (14) 



where rn H g = 2m h + k 2 jm^. Formal convergence of the integral is ensured 
by (slow) decrease of the Vg matrix element due to diminishing overlap be- 
tween the quarkonium bound-state wave function and the oscillating wave 
function of two heavy mesons in the continuous spectrum. However, the in- 
tegration over k which may be translated in the integration over the mass of 
the HH pair is to be cut-off at much smaller mg H values. Indeed, the consid- 
ered mechanism produces mass splitting of the 3 Pj states because of the mass 
difference between the pseudoscalar if-meson and its vector counterpart. We 
incur then that the scale of integration range is set by the m#* — mg mass 
difference. 

Results of numerical calculations confirm the qualitative conclusions drawn 
above. To investigate sensitivity of the effect we cut off the integration in (|TjJ) 
at some //(//> 2m//). Dependence of the .R-ratio for the 2P and IP states of 
the bottomonium and for the IP states of the charmonium are displayed in 
fig. 1 as a function of the reduced mass band r = (fJ>/2 — run) /(wis*— tuh)- As 
the initial values for i?'s we used those calculated in a potential model |]10| 



(using other predictions of potential models does not change qualitatively 
the obtained results). Note the fast onset of the R{2P) < R(1P) regime: 
already for r > 0.2 the hierarchy of the R values corresponds to the exper- 
imental observations (0) and it maintains with further increase of r. The 
considered mechanism does not yield quantitative agreement with the ex- 
perimental data: the overall set of data is shifted downwards with respect 
to experimental values (shown by error bars). A possible explanation may 
be that parameters of the spin-dependent potential determining the "initial" 
values of R were varied already in an attempt to reproduce experimentally 
observed regularities. It looks plausible that the consistent description of 
spin-dependent mass splitting in the quarkonium systems requires the ac- 
count of both the spin-dependent part of the Q — Q potential and of the 
bound-state coupling to mesonic channels. The investigation of the problem 
following these lines is in progress. 
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Summarizing, the process of (virtual) light-quark pair production con- 
siderably affects masses of the quarkonium states. Account of the emerging 
mass splitting of the 3 Pj triplet allows to reproduce qualitative features of 
the experimentally observed quarkonium level spacing. 

The author is indebted to F.Lev and S.Romanov for helpful comments. 
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